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ABSTRACT 


Models  for  deagglomeration  of  aggregates  of  particles  by  compressible 
flow  systems  and  for  the  fracture  of  particles  by  impact  grinding  have 
been  developed,  leading  to  a  prediction  of  particle-size  distributions, 
i.e.,  quasi -logarithmico-normal  distributions,  with  a  mean  geometric 
standard  deviation  of  1.6  ±  0.1. 

For  aerodynamic  deagglomeration,  the  size  distribution  is  related  to 
the  compressible  flow  properties  and  intrinsic  properties  of  the  particu¬ 
late  solids.  Impact  grinding  is  shown  to  be  governed  by  impact  probabil¬ 
ities  proportional  to  the  square  of  the  Stokes'  equivalent  spherical 
particle  diameter.  Data  are  presented  in  support  of  the  models,  and 
departures  from  the  models  are  discussed. 
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I .  INTRODUCTION 


The  aerosol  generation  of  particulate  solids  consists  either  of  fractur¬ 
ing  larger  particles  and  creating  an  air  dispersion  of  their  fragments,  or 
of  separating  preformed  particles  held  together  in  aggregates  by  surface 
contact  or  other  forces.  In  the  first  case,  new  surface  area  is  created, 
requiring  in  general  larger  expenditures  of  energy  and  applications  of 
larger  stresses.  Concomitant  with  the  exposure  of  fresh  surface  to  the 
gaseous  environment,  opportunity  for  catalysed  chemical  reactions  is 
afforded,  which  could  shift  the  powder  into  a  non-equilibrium  condition. 
Because  of  the  relatively  high  energy  and  stress  requirements  for  increas¬ 
ing  surface  area,  aerosols  of  powders  are  usually  generated  from  a  pre¬ 
sized  product,  prepared  either  by  grinding  or  spray  drying.  Special 
techniques  do  exist,  however,  for  the  preparation  of  certain  dry  materials 
that  make  possible  the  stable  aerosol  generation  of  non-presized  powders, 
through  reduction  of  stress  required  in  the  creation  of  new  surface  by 
alteration  of  powder  properties. 

In  general;  however,  any  particulate  solid  aerosol  generation  process 
involves  both  the  creation  of  fresh  surface,  or  grinding .  and  the  separation 
of  basic  presized  particles  held  together  as  aggregates,  or  deagglomeration . 
The  relative  extent  of  grinding  and  deagglomeration  depends  upon  the  energy 
available  in  the  generating  system  as  well  as  its  compressible  flow  charac¬ 
teristics,  the  particle-size  distribution  of  the  pre-dispersed  powder,  the 
design  of  the  aerosol  generating  mechanism,  and  the  intrinsic  .physical 
properties  of  the  solid  medium  comprising  a  particle. 

The  purpose  of  this  paper  will  be  to  present  models  for  grinding,  and 
for  particulate  solid  aerosol  generation  involving  compressible  fluid 
dynamics,  in  an  attempt  to  correlate  properties  of  the  generating  systems 
and  the  pre-dispersed  powder,  hopefully  leading  to  a  prediction  of  the 
initially  generated  aerosol  size  distribution.  Modification  of  the  initial 
aerosol  distribution  by  coagulation,  preferential  sedimentation  or  diffu¬ 
sion,  etc.,  will  not  be  discussed. 
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II.  AERODYNAMIC  DEAGGLOMERATION  OF  PARTICULATE  SOLIDS 


A.  CONFORMATION  OF  AN  AGGREGATE 

In  the  discussion  which  follows,  it  will  be  assumed  that  the  basic 
particles  comprising  the  powder  are  homogeneous  in  size,  a  basic  particle 
being  defined  as  the  smallest  volumetric  subdivision  of  the  powder  possible 
without  the  creation  of  fresh  surface.  This  assumption  does  not  strictly 
apply  to  most  powders  dealt  with  experimentally,  but  if  geometric  standard 
deviations  of  the  distributions  are  small  the  assumption  does  not  introduce 
any  serious  error.  A  proper  statistic,  such  as  Sauter  Mean  diameter,  aver¬ 
age  volume  diameter,  area  mean  diameter,  etc.,  can  be  chosen  to  approximate 
an  equivalent  "homogeneous  size"  for  a  given  application. 

It  will  be  further  assumed  that  the  basic  particles  are  initially 
separated  from  each  other,  i.e.,  no  aggregates  are  initially  present,  and 
that  they  are  in  random  motion  with  respect  to  each  such  as  applies  to  a 
turbulent  gas  suspension  or  to  a  mechanically  stirred  array.  Random  motion 
will  produce  collisions,  collisions  will  produce  aggregates,  and  a  steady 
state  will  be  reached  when  deagglomerating  forces  brought  to  bear  on  an 
aggregate  are  equal  to  forces  of  cohesion  holding  the  aggregate  together. 
Aggregates  consisting  of  various  numbers  of  basic  particles  can  be  formed 
in  two  limiting  ways: 

1 .  Type  1 

Type  1  results  from  random  collisions  between  basic  particles 
simultaneously  in  every  volume  element  of  the  array,  and  continuing 
collisions  between  resulting  aggregates,  etc.  Non-uniform  packing  of 
basic  particles  results. 

2.  Type  2 

Type  2  is  caused  by  the  initial  formation  of  an  agglomerate  by 
collision  of  two  basic  particles,  and  successive  build-up  of  this 
agglomerate  by  successive  collisions  between  it  and  individual  basic 
particles.  Uniform  packing  of  basic  particles  results. 

A  schematic  representation  of  agglomerate  formation  in  both  models, 
assuming  collisions  between  two  particles  for  illustration,  is  given  in 
Figure  1. 
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Figure  1.  Limiting  Types  of  Agglomerate  Formation. 
(Two-Particle  Collision  Assumed.) 
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If  Pc  Is  the  probability  of  a  single  1  collision  involving  two  particles, 
taken  as  equal  for  both  models,  then  the  probability  of  forming  an  agglom- * 
erate  consi&ting  of  basic  particles  for  each  type  is: 


I  lo«io  Ni 

p  m  f  6.301 
(N1>i  c  ' 

where  ■  2^ 

i  =  array  index 


(1) 


(2) 


The  ratio  of  probabilities  of  forming  an  agglomerate  of  type  1  to  type  2 
consisting  of  the  same  number  of  particles  is  then 


R  - 


g(«i) 1 
f<»l>2 


log  Nj. 
0.301 


+ 


(3) 


Since  f  <  1,  it  is  seen  that  aggregates  formed  by  the  mechanism  of 
type  1  are  much  more  probable  than  those  of  type  2.  For  example  for  Pc  *  0.5 
and  ■  8,  then  R  ■  16;  or  Pc  *  0.2,  ■  8  then  R  ■  625.  As  Pc  decreases, 

or  increases,  R  increases. 

Hence,  the  conformation  of  what  might  be  considered  a  most  probable  or 
representative  agglomerate  will  be  that  determined  by  probabilities  associ¬ 
ated  with  type  1,  having: 

(1)  a  quasi -spherical  shape,  due  to  the  spherical  symmetry  of  random 
motion  of  colliding  particles. 


(2)  five  sub-agglomerate  components,  on  the  average  being  the  smallest 
number  of  colliding  particles  defining  the  spherical  symmetry.  Each  of 
these  component  sub-agglomerates  will  consist  of  five  smaller  sub-agglomer¬ 
ates,  etc.,  ending  with  a  sub -agglomerate  consisting  of  five  basic  particles, 
which  would  correspond  to  the  first  aggregate  formed  in  the  initially 
dispersed  random  array. 


A  schematic  diagram  of  this  hypothetical  characteristic  aggregate  is 
shown  in  Figure  2. 
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The  largest  characteristic  aggregate  encountered  in  a  bed  of  pre¬ 
dispersed  powder  will  be  that  determined  by  the  magnitude  and  efficiency 
of  application  of  the  disrupting  stress  during  the  powder's  history, 
schematically  represented  in  Figure  2*  The  process  of  aerosol  genera¬ 
tion  then  can  be  considered  as  a  sequential  process  of  breaking  an  agglom¬ 
erate  of  type  1  into  its  component  sub -agglomerates,  etc.,  until  either 
stable  sub-agglomerate  or  basic  particles  are  dispersed.  An  important 
property  of  the  sub-agglomerate  hypothesis  is  the  decrease  in  numbers  of 
contacts  per  unit  area  between  component  particles  as  the  agglomerate 
increases  in  size,  i.e.,  as  the  agglomerate  gets  larger,  its  shear  strength 
or  tensile  strength  gets  smalleri  This  would  not  be  the  case  for  agglomer¬ 
ates  of  type  2  formation. 


B.  DISRUPTION  OF  AGGLOMERATES  BY  AERODYNAMIC  DRAG 

In  general,  mechanical  forces  capable  of  deagglomeratj.ng  particles 
arise  in  seven  ways  in  a  compressible  flow  system: 

1.  Impact  of  particles  with  fluid  flow  boundary  walls  due  to  vorticity 
of  the  particles'  motion.  This  could  result  from  high-order  turbulence  in 
the  fluid  flow  system,  dr  centrifugal  force  applied  to  a  particle  as  a 
result  of  its  motion  in  a  laminar  vortex.  The  average  velocity  of  the 
particle  in  this  system  would  be  in  general  parallel  to  the  impact  surface. 

2.  Impact  of  particles  with  each  other. 

3.  Centrifugal  disruption  of  a  particle  due  to  rotation  about  its 
own  axis. 

4.  Direct  impact  with  a  fixed  bounding  surface,  not  requiring  vorticity, 
where  the  average  velocity  vector  of  the  particle  intersects  the  impact 
surface. 


5.  Friction  drag,  or  shear,  resulting  from  velocity  gradients  in  a 
viscid  fluid  system;  and  form  drag,  or  pressure. 

6.  Inertia,  causing  preferential  acceleration  by  drag  forces  acting 
unevenly  over  particles  in  contact  and  composing  an  aggregate. 

7.  Explosive  decompression  of  an  agglomerate  containing  trapped  air 
in  its  voids  if  surface  pressure  is  rapidly  reduced.  This  effect  can 
theoretically  exist  for  a  spherical  agglomerate  under  drag  in  a  free 
stream. 

For  the  interaction  of  a  free,  unbounded  gas  stream  and  a  particle, 
mechanisms  (5),  (6)  and  (7)  would  predominate.  Drag  includes  both  a 
pressure  stress  and  a  shear  stress,  and  for  a  sphere  in  steady  laminar 
flow  of  a  viscuous  incompressible  medium,  or  for  a  compressible  medium 


at  Mach<il3,  these  stresses  are  normal  to  each  other  and  normal,  and 
parallel  to  the  surface  respectively.  The  resultant  drag  force  under 
these  conditions  is  parallel  to  the  direction  of  motion,  and  is  given  by 
Stokes'  law  for  spheres  with  radii  >  lp..  [The  Cunningham  correction  factor 
to  account, for  slip  is  required  for  radii  <  lp.] 

Ds  -  6  *  7?  .U  r  (4) 

U  *  relative  stream-particle  velocity 
r  «  radius  of  sphere,  or  characteristic  dimension 
=  viscosity  of  medium 

Stokes*  equation  holds  for  laminar  flow,  which,  for  a  sphere,  requires  a 
Reynolds  number,  Re  <  1 

Re  “  <  1 .  (5) 

p  ■>  density  of  medium 

A  general  expression  for  the  drag  or  resultant  force  on  a  particle  of 
any  shape  in  either  incompressible  or  compressible  media  for  either 
turbulent  or  laminar  flow  has  been  formulated  in  terms  of  a  drag  coeffi¬ 
cient.  Drag  acts  on  a  particle  in  direction  of  relative  motion  of  the 
medium  with  respect  to  the  particle. 


U  *  free  stream  velocity 


S  *»  projected  area  of  body  normal  to  stream 
Cp  ■  drag  coefficient 
p  =  density  of  medium  at  particle  surface 

The  drag  coefficient  is  a  function  of  shape  and  Reynolds  number,  and 
for  irregular  shapes  is  empirically  determined.3  A  function  relating  Cp 
and  Rg  for  spheres  over  the  range  of  0  <  Rg  <  2000  has  been  derived  by 
Langmuir3  with  average  deviations  from  experimentally  determined  values 
of  less  than  4  per  cent. 
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CD  Re/24  -  1  +  0.197  Re°-63  +  2.6  x  10~4  Rg1,38 


(6)A 


If  a  particle  is  free  to  conform  to  an  external  pressure  stress  over 
its  surface,  its  shape  will  be  altered  if  the  resultant  stress  is  not 
uniform  in  magnitude  and  direction  over  every  region  of  the  surface.  It 
can  be  shown  that  particles  obeying  Stokes'  law  are  not  distorted.  How¬ 
ever,  particles  accelerated  in  compressible  media  do  experience  non-uniform 
as  well  as  non-symmetrical  pressure  stresses,  and  aggregates  of  particulate 
solids  or  liquid  droplets  will  deform  during  acceleration  if  the  magnitude 
of  the  stress  exceeds  the  tensile  stress  or  surface  tension  pressure  of  the 
respective  particles.  • 

For  a  sphere  under  acceleration  in  an  inviscid  fluid  at  rest  at 
infinity,  the  surface  pressure  distribution  as  derived  from  a  velocity 
potential  can  be  shown  to  be4 

F-~-£p.  -  I  r  Cos  6  +  4  U2  (9  Cos20  -5)  (7) 

p  2  8 

where  r  *  radius  of  sphere 

0  *  angle  measured  positively  from  flow  axis  from  stagnation  point 

U  *  velocity  of  sphere  in  inertial  frame  of  reference 

P  =  pressure  on  surface 

PQ  =  stagnation  pressure 

t  =  time 
2 

Since  the  Cos  0  term  is  symmetrical  about  the  axis  of  motion,  the  pressure 
distribution  contributing  to  drag  is : 


P  ~  Po  =  i  r  Cos  0 
P 


dU 

3t 


(8) 


Equation  8  is  plotted  in  Figure  3  for  the  case  r  *  12|i  and  an  initial  Mach 
0.1  relative  particle  free  stream  velocity  at  atmospheric  pressure.  At 
Mach  0.1,  compressibility  effects  are  negligible.  From  Newton's  second 
law  and  Equation  (6) 
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particle  density 


The  radial  pressure  distribution  of  Equation  8  is  asymmetrical  and 
non-uniform,  thereby  contributing  to  acceleration  of  the  particle  in  an 
inertial  frame  and,  to  some  extent,  to  flattening  of  a  particle  in  the 
direction  of  relative  fluid-particle  motion  due  to  unequal  acceleration 
of  its  fore  and  aft  surfaces. 


The  second  term  in  Equation  (7),  however,  for  the  conditions  of  flow 
defined  above,  is  greater  than  Equation  (8)  by  orders  of  magnitude,  and, 
although  not  contributing  to  drag  because  of  its  symmetry,  does  have  the 
effect-  of  producing  greater  particle  distortion  in  a  body-fixed  frame  of 
reference.  This  "deforming"  pressure  distribution  is  plotted  in  Figure  4. 
It  is  seen  that  the  minimum  pressure  is  sub-atmospheric  at  0  >  90°. 

A  rough  estimate  of  the  rate  of  deformation  of  flattening  for  the 
case  under  consideration  can  be  made  by  the  following  considerations: 
Consider  a  unit  area  on  the  surface  at  0.  Then  acceleration  of  the  surface 
of  an  incompressible  inviscid  sphere  in  a  body-fixed  frame  is  given  by 

p  U2  2  d2x 

(P  -  PQ)  =  (9  Cos20  -5)  -  Ps  ^7  (10) 


Ps  *  surface  density  of  particle 
X  =  radial  distance  moved  by  surface  element,  0  constant 

giving 

P  U2  2  2 

x0 =  rr (9  cos  9  _5)  t2  (ll) 

but 


X0  "  ro  *  r 

rQ  =  initial  spherical  radius 
r  -  distorted  radius 


Hence,  for  0  =  Cf  and  Xg  =  rQ,  then  t  *  4.5  x  10"^  sec.  Or,  for  0  =  90° 
and  t  *  4.5  x  10 sec.,  then  i-g  m  ggo  ■  j-  rQ  or  the  sphere  tends  to 
develop  into  a  crenated  disc  with  a  1257.  increase  in  radius.  The  pressure 
distribution  will  approach  that  over  a  disc  with  its  plane  normal  to  flow. 
This  pressure  distribution  over  the  leading  surface  is  given  by6 


P  =  P  - 
rr  co 


a2  r2 


2 


(12) 


RELATIVE  FLUID 
VELOCITY 


r 


Figure  4.  Initial  Deforming  Pressure  Distribution. 


where 


r  =  radius  vector 


Pq  =  stagnation  pressure 
a  =  constant 

Cavitation  will  tend  to  develop  behind  the  sphere  as  it  distorts,  producing 
an  asymmetrical  pressure  distribution,  and  thus  will  further  increase  the 
drag.  Hence,  when  deformation  of  the  sphere  begins,  disruption  will  result 
if  the  applied  stress  is  maintained. 


The  aerodynamic  pressure  stresses  and  viscous  shear  stress  would  then 
probably  play  the  following  sequential  roles  in  the  deagglomerating  process 
of  a  spherical  aggregate: 

(1)  Deformation  or  flattening  of  the  sphere  by  normal  components  of 
stresses,  Figure  5a. 

(2)  Further  deformation  due  to  normal  aerodynamic  pressure  stress  on 
the  resulting  concave  surface,  causing  continued  radial  separation  of  the 
aggregate  with  respect  to  its  center  of  mass.  As  the  sphere  flattens,  a 
highly  non-symmetrical  pressure  distribution  occurs  that  will  further 
accelerate  the  particle.  If  the  pressure  stress,  greatest  at  the  stagna¬ 
tion  point,  exceeds  the  tensile  strength  of  the  aggregate,  it  will  begin 
tearing  through  the  axis  of  symmetry,  Figure  5b. 


(3)  The  tangential  shear  stress  due  to  gradients  in  the  boundary  layer 


will  produce  an  additional  radial  component  of  acceleration  of  the  surface, 
Figure  5c.  For  a  lamina^  boundary  layer  (Re  <  3.000)  over  a  disc,  this 
shear  stress  is  given  by 


(4)  The  inertial  force  of  the  aggregate  acts  through  its  center  of 
mass  in  a  body  fixed  frame,  Figure  5d. 

Hence,  during  acceleration  of  the  aggregate,  a  couple  is  produced  with 
radial  components.  The  net  effect  is  to  increase  the  axial  stress  parallel 
to  flow,  with  radial  acceleration  away  from  the  center  of  mass.  The  shear 
stress  on  the  surface  produces  additional  relative  motion  of  the  surface 
with  respect  to  the  interior  of  the  aggregate,  which  moves  with  a  radial 
component  because  of  the  off-axis  normal  pressure  stress.  Deagglomeration 
results  from  this  relative  motion  of  the  component  parts,  or  sub-aggregates, 
away  from  each  other. 

Another  possible  factor  contributing  to  deagglomeration  is  "cavitation" 
in  the  agglomerate  voids  when  the  static  pressure  on  the  surface  falls  below 
the  ambient  pressure  inside.  A  maximum  pressure  gradient  for  the  potential 
flow  case  considered  would  occur  at  0  *  90°,  leading  to  an  "explosive  decom¬ 
pression"  radially  and  normal  to  the  relative  stream  velocity. 


20 


III.  THEORETICAL  AEROSOL  SIZE  DISTRIBUTIONS 


A.  MATHEMATICAL  FORM 

1.  Drag-Induced  Deagglomeration 

It  was  shown  above  that  an  agglomerate  consisting  of  small  uniformly 
sized  basic  particles  with  non-uniform  packing  is  more  probable  than  one 
with  uniform  packing.  Further,  from  symmetry  consideration  of  random  colli¬ 
sions,  a  quasi-spherically  shaped  agglomerate  consisting  of  5  sub-agglomer¬ 
ates  would  be  the  most  probable  non-uniform  array.  Hence,  an  aggregate  is 
considered  to  consist  of  sub-agglomerates,  each  in  turn  of  sub-agglomerates, 
etc.,  and  the  deagglomeration  process  to  consist  of  breakup  of  a  particle 
into  its  composing  sub-agglomerates  in  a  single-cycle,  step-wise  process. 

An  aerosol  generated  from  a  mass  of  uniformly  sized  basic  particles 
will  have  an  upper  bound  on  its  Bize  distribution,  i.e.,  a  largest  size 
will  exist,  and  all  aggregates  larger  than  this  size  will  have  been  deagglom- 
erated.  Deformation  of  an  agglomerate  will  probably  take  place  when  an 
applied  aerodynamic  "deforming"  stress  exceeds  the  minimum  tensile  stress 
of  the  aggregate,  and  once  deformation  occurs,  deagglomeration  follows.  A 
measure  of  the  probability  of  deformation  and  therefore  of  deagglomeration 
of  an  aggregate  into  its  sub -aggregates  can  be  expressed  by: 


n 


P  *  probability  of  deagglomeration 

D  3  available  deagglomerating  force,  assumed  constant  for  a 
given  process 

Fc  =  minimum  cohesive  force  holding  aggregate  together 

K  =  factor  of  proportionality,  constant  for  a  given  process 

For  aerodynamic  drag, 

2  2 
D  =  Cp  P  U  n  d4 

8 


(14) 


(6) 


d  =  diameter  of  spherical  agglomerate 
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Also 

Fc  =  0fc  (15) 

P  =  number  of- contacts  between  sub -agglomerates  in  plane 
through  center  of' agglomerate 

fc  -  force  per  contact  between  sub-agglomerates 

But 

p  =  (n-2)3  (16) 

n  «=  number  of  sub -agglomerates  in  each  agglomerate  (n  =  5  as  a 
probable  value) 


And  for  osculating  spherical  svfcffigglomerates  where  point  contact  is  assumed 


fcj,  *  force  per  contact  between  two  basic  particles 
Therefore 


Hence 


Fc  “  <n  "  2>3  fcb 

K  Cn  P  U2 
P  =  8(n-2)3  fcb 
nd? 


(17) 


(18) 


(19) 


For  all  particles  equal  to  or  larger  than  the  upper  bound  on  the  aerosol 
size  distribution,  which  have  been  broken  down  by  the  application  of  the 
deagglomerating  stress 


P 


K  Cp  P  U2 

=  3(n-2)3  =  1  (20) 


d  =  upper  bound,  or  smallest  sized  aggregate  with  a  probability 
of  deagglomeration  of  unity 

And  for  a  given  deagglomerating  process,  characterized  by  a  constant  drag 
force 


8  (n-2)3  fcb 
*  Cp  P  U2d^2 


(21) 
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K  can  be  taken  as  a  relative  measure  of  the  efficiency' with  Which  available 
specific  energy  in  the  gas  stream  is  utilized  in  effecting  application  of 
the  deagglomeration  stress. 

Combining  (19)  and  (21),  the  probability  of  deagglomerating  any 
other  aggregate  of  diameter  d^  <  d^,  into  its  composing  sub -aggregates, 
is  given  by 


P.  = 

X 


(22) 


But 


d-£  =  diameter  of  aerosol  agglomerates 
i  =  index  denoting  stage  of  successive  agglomerate  break-down, 
*’  0,  1,  2  — 


For 


-  d0  +  e 


(23) 


where  dQ  =  representative  number  weighted  diameter  of  largest  particle 
in  aerosol,  after  deagglomeration  process  is  complete 

e  =  arbitrarily  small  number 


e  -  0.005  d. 


Pi  = 


0.99  dt2 


(24) 


But  for  n  =  5 

“i+i  •  f 


(25) 


There  fore 


,  0.99 

Pi  =  22i 


(26) 


The  number-size  distribution  of  a  stable  aerosol,  after  the  dissemi¬ 
nation  process  is  complete,  can  be  computed  on  the  basis  of  Equation  (26), 
(25),  and  (23).  A  series  expression  with  histogram  class  intervals,  rather 
than  a  continuous  distribution  function,  results.  Random  fluctuations 
about  mean  values  for  each  class  interval  would  physically  result  in  a 
continuous  distribution,  however. 
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The  number  of  particles  in  each  size  fraction  of  the  generated 
aerosol  is  computed  from  the  difference  between  the  number  entering  a 
size  class  (from  the  size  class  immediately  preceding  it  in  the  direction 
of  increasing  size)  and  the  number  leaving  it  as  a  result  of  the  next 
successive  stage  in  the  deagglomeration  process. 


7f0  =  No  (1-PQ)  — .  (27) 

i  =  0 

JJ  =  number  particles  in  zerot^1  class  corresponding  to  size  dQ 
No  =  total  number  of  particles  entering  class  with  d  =  dQ 
Similarly 

7?!  =  <NonVNo  n  P'o  Pl>  <28> 

i  =  1 

n  =  representative  number  of  sub-agglomerates  in  an  aggregate 

n  constitutes  the  average  number  of  pieces  into  which  a  particle  will  break 
as  a  result  of  a  single  stage  in  the  deagglomeration  process,  taken  as 
5  in  the  agglomerate  model.  Hence,  the  general  expression  for  the  number 
of  particles  in  the  generated  aerosol,  as  a  function  of  the  deagglomeration 
stage,  is  given  by 

\  ‘  N0  n1  Pi-l  Pi-2-— Po  <l  ’  Pi>  <29> 


or,  using  Equation  (26)  for  P^, 

7Ji  «  N0  P0  2-1  <zi  -  Po  2_i> 


(30) 


The  average  diameter  of  an  agglomerate  associated  with  each  class,  i,  in 
Equation  (30)  is  obtained  from  Equation  (26). 


d 


i 


(31) 


i  *  0,  1,  2  — 


dD  is  determined  by  the  physical  parameters  of  the  deagglomeration  process 
from  Equations  (16),  (20),  (23) 


do 


/  £  P  ^cb 

Ik  j:  c_.p  u2 
\  - D 


% 


(32) 
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A  continuous  cumulative  number-size  distribution  can  be  constructed  by 
summing  over  i  and  associating  7].  with  sizes  of  particles  having  an  upper 
class  value  of  d, .  Number-size  distributions  for  the  case  of  n  =  5  are 
shown  in  Figure  6  with  various  values  of  dQ  as  a  parameter.  A  quasi -log 
normal  distribution  results  with  an  average  standard  deviation  of  1.6  ±  0.1. 
The  slope  is  invariant  with  the  value  assumed  for  n. 

The  distribution  is  of  course  bounded  by  dQ  and  d  ,  where  dg  = 
diameter  of  the  basic  particle.  This  precludes  an  explicit  log-normal 
distribution  form  as  mathematically  inferred  by  Kolmogoroff  for  particles 
resulting  from  grinding.  It  differs  significantly  from  that  of  Martin 
and  subsequently  derived  by  Dallavalle,  who  found  the  number-size  distri¬ 
bution  function  in  the  case  of  fracture  to  follow  the  compound  interest  law.* 


The  transformation  from  a  cumulative  number-size  distribution  to  a 
cumulative  weight  distribution  can  be  accomplished  in  two  ways.  Constant 
particle  density  is  assumed. 

1.  Assume  a  quasi-log  normal  fit,  and  use  an  average  value  of  a 
for  the  two  straight-line  portions  of  the  curve: 

Then 

d  =  d  e  d  ■  number  median  diameter  (33) 

cm  c 


CTg  =  1.60  ±  0.10 


(34) 


or 

d  ®  1.95  d'A 
m  c 


d  »  mass  median  diameter 
m 


(35) 


2.  Make  no  assumptions  concerning  the  theoretical  form  of  the 
distribution,  and  compute  the  distribution  by 


f^m>i  "  f(n)i 


(36) 


The  increase  in  slope  of  the  cumulative  number-size  distribution 
between  -£■  and  dQ  is  due  to  the  averaging  effect  of  the  model,  which 
precludes  non-integral  changes  in  agglomerate  sizes.  For  this  large  size 
class  the  series  summation  does  not  adequately  reflect  a  continuum.  Con¬ 
sequently,  method  number  1  for  transforming  the  number  distribution  to  a 
mass-weighted  distribution  is  the  more  accurate.  Variations  of  agglomerate 
density  with  size  have  been  neglected  in  the  transformation.  d0  corresponds 
to  the  97%  in  the  mass  transformed  distribution. 


*  Georgia  Institute  of  Technology,  Atlanta,  Georgia.  "Studies  and  investi¬ 
gations  of  agglomeration  and  deagglomeration  of  solid  particles,"  by 
C.  Orr  and  J.M.  Dallavalle.  Semifinal  Report.  30  June  1956.  Contract 
DA- 18-065 -404-CML-88 . 
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Theoretical  Cumulative  Number-Size  Distributions  of  Dry  Aerosols 
With  Maximum  Representative  Size  •  do>  as  Parameter. 
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Hence,  in  summary: 

For  number-size  distributions 
dc  =  0.21  dQ 
ag  =  1.60  ±  0.1 

dQ  =  corresponding  to  the  99.99  percentile  on  a  logarithmic 
probability  plot 

For  mass-size  distributions 

d  =  0.42  d„ 
in  o 

CTg  =  1.60  ±  0.01 

dQ  =  corresponding  to  the  97.0  percentile 

A  plot  relating  the  median  size  and  the  size  corresponding  to  the  15 
percentile  to  the  representative  largest  size  present  in  the  aerosol,  dQ, 
is  given  in  Figure  7  for  both  number-  and  mass-weighted  distributions. 

The  vertical  lines  on  the  mass  distribution  reflect  the  spread  of  values 
between  the  two  methods  of  computing  the  transformation.  The  curve 
corresponds  to  method  1. 

2.  Size  Distribution  of  Particles  by  Impact  Grinding 

The  size  reduction  of  basic  particles  by  impaction  grinding  under 
conditions  of  turbulent  flow  can  be  considered  to  consist  of  two  independ¬ 
ent  mechanisms,  both  associated  with  a  probability,  and  the  probability  of 
breakup  can  be  computed  on  the  basis  of  the  product  of  probabilities  of 
each  mechanism  occurring.  The  mechanisms  are:  impact  with  a  surface  as 
a  result  of  a  particle  leaving  a  stream  line,  and  fracture  as  a  result  of 
the  impact. 

Hence 

-*l'f  (37) 

P  =  probability  of  breakup 
B 


p 

i 


-  probability  of  impact 


P^  -  probability  of  fracture 


Diameter,  u,  at  Stated  Percentile 
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The  probability  of  a  particle  leaving  a  curved  stream  line  and  colliding 
with  a  bounding  surface  is  proportional  to  the  centrifugal  force  acting 
on  the  particle  due  to  its  inertia  and  the  shear  drag  opposing  radial 
motion 


"5” P  d3  U2 
P,  =  K  - E - 


3  «  77  0  d 


(38) 


or 


Pi  = 


i  2 

K  P^U  dZ 


1  18^  R 


(39) 


R  =  average  radius  of  curvature  of  stream  lines  at  impact 
sur  face 


K  =  constant  of  proportionality 


The  probability  of  fracture  upon  impact  can  be  considered  proportional  to 
the  specific  energy  available  for  fracture  work  and  inversely  proportional 
to  the  specific  energy  required  for  fracture.  Particles  undergoing  plastic 
deformation  would  tend  not  to  fracture. 


P 


f 


i  t 


K 


*Ppu2 

y  * - 

S  /2E 


(40) 


S  =  critical  compressive  stress  of  particle 
E  *  modulus  of  elasticity  of  particle 
Hence,  probability  of  breakup  under  impact,  using  Equation  (37) 

|  i  9  3  5  9 

K  P_  U  <r  K  PD  UZE 
p_  = - E - _E - 


B 


(41) 


18  77  R  S 


Or  for  a  given  system 


(42) 


K1 ' '  =  constant 


Data  have  been  obtained  by  Epstein®  in  support  of  Equation  (42)  for  the 
impact -breakage  probability  being  proportional  to  d*,  in  contradiction  to 
a  basic  assumption  in  his  derivation  of  logrithmico-normal  distribution 
for  breakage  of  solids. 


? 
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Applying  the  same  reasoning  used  in  the  derivation  of  Equation  (30), 
where,  in  place  of  an  upper  characteristic  aggregate  size  an  initial  uniform 
basic  particle  size  exists,  an  identical  distribution  of  numbers  of  parti¬ 
cles  as  a  function  of  stage  in  the  breakup  process  is  obtained.  Sizes  are 
quasi-log  normally  distributed,  if  the  sizes  related  to  successive  stages 
are  in  a  constant  ratio. 

For 

di 

-  **  2 .  then  <7  =  1.6  ±  0.1 

di+i  8 

The  geometric  standard  deviation,  c_,  depends  on  the  ratio  of  sizes  but 
is  independent  of  n,  the  number  of  fragments  resulting  from  a  single 
fracture.  For  spheres,  n  ■  4,  or  5;  for  flat  platelets,  n  «  2,  if 
dt 

di+l  "  2‘ 


B.  EXPERIMENTAL  OBSERVATIONS 

1.  Aerodynamic  Deagglomeration 

A  major  difficulty  in  determining  aerosol  properties  is  the  effect 
of  the  measuring  device  on  the  parameter  of  interest.  This  is  particularly 
true  of  devices  that  collect  a  "representative  sample"  of  the  aerosol,  such 
as  the  Cascade  impactor10  for  example.  Light -scattering  techniques,  such 
as  those  of  Dimmick1  and  O'Konski,13  though  they  do  not  require  the  collec 
tion  of  representative  samples  and  hence  do  not  distucb  the  system  being 
measured,  nonetheless  require  careful  interpretation,  involving  a  knowledge 
of  the  interaction  of  the  containing  chamber  with  the  aerosols,  as  well  as 
Interpolation  of  light -scattering  phenomena.  A  further  complication  is 
the  inherent  instability  of  an  aerosol  due  to  collision  processes  involv¬ 
ing  the  suspended  particles,  i.e.,  random  collisions  due  to  turbulence 
or  Brownian  motion.13 


One  system  that  tends  to  minimize  sampling  and  interpretation 
problems,  and  that  is  suited  for  aerodynamic  aerosol^ generation  studies, 
has  been  found  to  be  the  Sharpies  Micromerograph.14’  Designed  for 
small  quantities  of  material  (<  50  mg)  it  measures  a  cumulative  mass-**- 
Stokes'  drag  diameter  distribution  of  the  entire  aerosol,  and  has  been 
found  satisfactory,  with  careful  handling,  for  the  determination  of 
particle  sizes  >  1  p  d.  It  operates  on  the  principle  of  aerodynamic 
drag,  subjecting  a  smair~mass~dT  powder  to  a  steep-fronted  sonic  air 
blast  at  controllable  stagnation  pressures.  The  particulate,  solids 
entrained  in  the  sonic  flow  pass  through  an  expanding  conical  orifice 
or  nozzles  of  special  design,  and  the  steady-state  aerosol  is  measured 
by  a  continuous  cumulative  weighing  of  particles  that  have  settled 
according  to  Stokes'  Law,  Figure  8. 


3 


3 


Figure  8.  Schematic  diagram  of  the  Hicromerograph 
Deaggloxaeratzion  Section. 
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One  series  of  studies  on  the  effect  of  gas  density  on  breakup* 
with  the  Micromerograph  used  as  the  test  powder  Serratia  marcescens,,  a 
ground,  freeze-dried  vegetative  microorganism  contained  in  an  inert 
crystalline-like  matrix.  It  had  a  mass  median  diameter,  dm,  ■  4.1  |i 
and  a  geometric  standard  deviation  of  =  1.6. 

Similar  studies,  using  capillary  tubes  as  dispersing  nozzles, 
were  conducted  by  Orr  and  Dallavalle  with  3.3  p.  MMD  S_.  marcescens 
product.  The  basic  size  distributions,  though  not  conforming  to  that 
strictly  required  of  the  model,  t.e.,  homogeneity  of  size,  small  median 
diameters  and  geometric  standard  deviations.  The  characteristic  basic 
particle  size  for  the  heterogeneous  distributions  would  correspond  to 
the  average  volume  diameter.  The  mass -size  distribution  data  obtained 
from  both  these  studies,  for  a  range  of  stagnation  pressures  of  50  psig 
to  400  psig,  have  been  compared  with  the  aerodynamic  aerosol  generation 
model  presented  above,  in  an  attempt  to  correlate  compressible  flow  and 
particulate  solid  properties  with  the  aerosol  size  distribution  as 
follows : 


On  the  basis  of  the  model,  the  absolute  size  distribution  is 
defined  in  terms  of  dQ,  the  representative  largest  agglomerate  present 
in  the  aerosol.  This  value  in  turn  is  related  to  the  particle  properties 
and  compressible  flow  properties  by  Equation  (32). 


d 


o 


I  8  P  fcb  \  % 

\  K  « s p  u2 ) 


(32) 


All  quantities  with  the  exception  of  K  are  either  measured,  calculated  or 
inferred  from  the  model. 

a.  dQ,  maximum  measured  representative  agglomerate  size, 
corresponding  to  the  97  percentile  on  a  mass  weighted  basis. 

b.  f  the  force  of  cohesion  between  sub -aggregates .  Assuming 
point  contact;,  this  is  equal  to  the  force  between  basic  particles.  A 
minimum  value  for  this  force  has  been  obtained  from  bulk  tensile  strength 
measurements  of  S_.  marcescens  powder,  assuming  close  packing  of  basic 
particles.**  Bulk  tensile  strengths  are  of  the  order  of  2000  dynes/cm^ 
for  a  packing  fraction  of <~'0.20. 

fck  >  1.8  x  10  ^  dyne 


*  Bures,  M.G.,  Derr,  J.S.,  unpublished  data. 

**  Derr,  J.S.;  Bures,  M.G. ;  Gordon,  G.,  unpublished  data. 


c.  U,  gas  velocity,  sonic  at  the  powder-gas  interaction  rate. 

d.  p,  gas  density,  computed  for  isotropic  flow  at  Mach  1,  for 
the  various  stagnation  pressures. 

e.  Cq,  drag  coefficient  computed  on  the  basis  of  Reynolds  number, 
with  d0  as  the  characteristic  dimension.  Values  ranged  from  0.38  to  0.50. 

f.  p,  number  of  contacts  between  sub-aggregates.  For  n  -  5 
p  =  9,  assuming  close  packing. 

Values  of  p/k  were  calculated  using  Equation  (32),  and  a  mean  was 
determined  for  all  experimental  conditions,  with  the  95  confidence  interval 
on  the  mean  expressed. 

p/k  =  2.2  x  104  (±  0.5  x  104  957.  C.I.)  (43) 

Using  this  value  of  p/k,  theoretical  mass-size  distributions  were 
computed  using  Equations  (30),  (31),  (32)  for  the  values  of  gas  density 
used  in  the  experiments.  A  mean  value  of  Cq  *  0.43  was  also  used,  and 
sonic  flow  assumed.  _ 

A  comparison  of  the  theoretical  distribution  with  experimental 
average  values  is  given  in  Table  I.  The  experimental  geometric  standard 
deviation  is  a  mean  of  the  50/16  and  84/50  percentiles,  as  is  the  theo¬ 
retical  GSD.  The  agreement  is  close,  considering  the  variety  of  conditions 
over  which  p/k  was  computed  and  the  departure  from  experimental  conditions 
of  the  model  assumption  of  homogeneity  of  basic  particle  size.  Generali¬ 
zation,  however,  cannot  be  inferred  from  the  specific  absolute  value  of 
p/k  obtained  from  these  experiments  without  additional  work  with  other 
powders.  As  the  mass  median  diameter  of  the  basic  heterogeneous  particle 
size  distribution  increases,  measured  values  of  p/k  for  a  given  fcb  would 
increase,  reflecting  a  larger  apparent  dQ  than  a  homogeneous  model  would 
predict.  What  is  perhaps  of  interest  is  the  apparent  constancy  of  the 

ratio  a  4  under  the  conditions  studied  for  the  S_.  marcescens  material. 

Iv 

Values  for  both  p  and  fc^,  are  minimal.  Therefore,  k"^  is  maximal. 

K"1  <  2.4  x  103  ±  (0.6  x  103) 

K  ^  is  analogous  to  the  critical  Weber  number,  Wec,  for  liquid 
droplet  instability  under  aerodynamic  drag.17’1* 

Wec  =  “  >  3  to  10  (44) 

D  =  drag  pressure 
r  =  drop  radius 


a  =  surface  tension 
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TABLE  I.  ABSOLUTE  MASS -SIZE  DISTRIBUTION  OF  S.  MARCESCENS 
POWDER -SONIC  DEAGGLOMERATION  MODEL 
COMPARED  WITH  EXPERIMENT 


Stagnation 

Pressure 

"4.1" 

S.  marcescens 

"3.3" 

S.  marcescens 

Theory 

50  psig 

c 

7.4  p 

9.0  p 

11.5  ±  3.2*  p. 

CTg 

1.61 

1.75 

1.60  ±  0.1 . 

100  psig 

c 

7.4 

8.0 

8.3  ±2.0. 

CTs 

1.68 

1.67 

1.60  ±  0.1 

200  psig 

dm+ 

6.7 

6.9 

6.1  ±  1.5 

CTg 

1.60 

1.71 

1.60  ±  0.1 

300  psig 

d  + 
m 

6.3 

5.8 

5.1  ±  1.4 

CTg 

1.58 

1.62 

1.60  ±  0.1 

400  psig 

<m+ 

5.8 

4.3 

4.5  ±  1.1 

°g 

1.51 

1.54 

1.60  ±  0.1 

*  95  per  cent  confidence  interval. 

+  Mass  median  diameter. 

Nozzles:  Conical  deagglomerator,  3  tubes  V  L  x  0.18",  0.08", 
0.04"  ID.  " 


The  apparently  high  ratio  of  magnitude  of  K  Vwe  may  well  be  due  to  the 
uncertainty  of  the  fcb  for  which  a  minimum  value  has  been  measured'.  The 
presence  of  30  p  d  aggregates  in  the  cleavage  plane  of  the  tensile  strength 
measurements  mentioned  above  would  allow  a  hundredfold  increase  in  fc^  and 
result  in  We  =  K"^.  Further  work,  however,  is  required  to  clarify  this 
comparison  with  liquid  droplet  breakup. 

Mass-size  distributions  determined  by  the  Mi cromer ograph  technique 
for  a  variety  of  other  powders  have  been  compared  with  the  shape  predicted 
by  the  model,  i.e.,  quasi-log  normal  distribution  with  a  geometric  standard 
deviation  of  1.6  ±  0.1.  These  are  shown  in  Figure  9.  For  those  powders 
with  basic  mass  median  diameters,  dm  <  2  p  as  determined  either  by  optical 
sizing  or  with  the  Whitby  Centrifuge  technique,19’ So  the  agreement  with 
prediction  is  good.  Sodium  fluorescein  and  sulphur  aerosols,  characterized 
by  dm  >  14  ft  and  basic  og  >  2,  depart  from  the  model.  Theoretical  curves 
were  fitted  to  the  experimental  ones  at  the  median  diameter. 


Magnesium  Silicate  50+ 


J  UIL 


Cumulative  Per  Cent  Mass  Thru  Stated  Size 
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For  those  powders  conforming  to  the  model,  fcb0/K  was  computed 
from  Equation  (32)  for  sonic  flow.  The  results  are  given  in  Table  II. 
Interpretation  of  these  data  must  await  further  work  on  measurements  of 
fcb«  It  is  interesting  to  note,  however,  that  there  is  order  of  magnitude 
agreement  for  the  various  types  of  solids,  and  in  the  case  of  magnesium 
silicate  close  agreement  in  the  values  despite  the  large  difference  in 
stagnation  pressure.  Since  fcb  and  (3  are  properties  of  the  powder,  and 
fcb(3/K  approximately  equal  for  the  two  pressures,  K"1  should  be  equal  and 
therefore  independent  of  gas  density.  The  same  conclusion  is  reached  in 
-the  absolute  size  distribution  correlations  discussed  above  for  the  S. 
marcescens  powder.  Hence,  it  may  not  be  unlikely  that  K-3  has  the  same 
significance  in  solids  aerosol  dissemination  that  the  critical  Weber 
number,  Wc,  has  for  liquid  aerodynamic  breakup. 


TABLE  II.  AERODYNAMIC  DEAGGLOMERATION  PARAMETER,  fcb0/K, 
FOR  VARIOUS  POWDERS -SONIC  FLOW 


Powder 

Stagnation 

Pressure 

(paig) 

Gas  Density 
(gm  cm-3  x  10”3) 

do 

(n) 

°D 

fcbP 

(dy&es) 

Magnesium  silicate 

50 

2.84 

31.0 

0.43 

4.9 

Zinc-cadmium-sulphide* 

200 

9.5 

24.0 

.46 

11.2 

Acrylic  resin  (FT)* 

400 

13.9 

15.0 

.43 

7.9 

Bacillus  subtilis* 

400 

13.9 

13.5 

CM 

• 

6.3 

Magnesium  silicate 

400 

13.9 

11.0 

.41 

4.1 

Cadmium  oxide* 

400 

13.9 

7.0 

0.41 

1.7 

*  Reference  17, 


2.  Fluid  Energy  Mill  Grinding 

Grinding  by  a  fluid  energy  mill  is  one  effective  way  of  reducing 
the  basic  particle  size  of  solids.  In  a  sense,  a  fluid  energy  mill  can 
be  considered  a  special  type  of  solids  aerosol  generator,  where  high  impact 
efficiency  and  shear  rates  are  made  possible  through  large  energy  inputs. 
The  larger  tensile  or  shear  strength  of  basic  particles  as  compared  with 
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basic  particulate  agglomerates  requires  increased  energy  and  efficiency  for 
grinding  as  compared  with  deagglomeration.  For  example,  for  £,.  marcescens 
powder  the  basic  solid  shear  strength  is  of  the  order  of  10?  dynes/cm^  as 
compared  with  2  x  10^  dynes/cm^  for  the  bulk  powder  tensile  strength  of 
4  |i  d,ji  product. 

By  controlling  fluid  energy  grinder  parameters,  it  is  possible  to 
reduce  basic  particle  size  in  an  empirically  controlled  way  for  a  given 
material.  Results  of  such  a  study  with  the  freeze-dried  S_.  marcescens 
powder  are  shown  in  Figure  10. sl  This  material,  characterized  by  crystal¬ 
line  structure  and  a  high  bulk  modulus  of  elasticity,  is  readily  grindable. 
The  basic  particle  mass-size  distributions  shown  in  the  figure  were  deter¬ 
mined  by  use  of  the  Whitby  centrifuge  technique,  representing  that  of  the 
completely  deagglomerated  powder.  Curve  A  resulted  from  an  initial  ball- 
milling  operation,  producing  particles  with  platelet  shape  having  a 
characteristic  thickness  of  »  5  (i.  Powders,  whose  distributions  are 
shown  in  curves  BC  and  D  of  Figure  10,  were  obtained  by  the  further  fluid 
energy  mill  grinding  of  this  ball-milled  product,  using  different  values 
of  the  grinder  parameters.  The  experimental  data  are  fitted  to  the  model 
predictions  at  the  median  value,  and  it  is  seen  that  logrithmico-normal 
mass-size  distributions  with  a  geometric  standard  deviation,  Og,  of  1.6  ± 
0.1  are  obtained,  as  inferred  by  the  discussion  in  Section  III,  A,  2. 


IV .  DISCUSSION 


Models  leading  to  the  prediction  of  initially  generated  particle  size 
distributions  of  particulate  solid  aerosols  and  ground  powders  have  been 
developed.  For  the  aerosol  case,  a  homogeneous  distribution  of  basic 
particles  is  assumed,  although  powders  with  small  median  diameters  and 
standard  deviations  afford  working  approximations.  Data  with  such  pmrders 
have  been  obtained  for  sonic  deagglomeration,  using  the  Mi cromer ograph 
under  a  variety  of  compressible  flow  conditions.  Under  the  conditions 
studied,  data  have  been  found  to  support  the  model  predictions  for  the 
shape  of  the  size  distribution,  i.e.,  quasi-log  normal  with  Og  of  1.6  ±  0.1. 

For  powders  having  median  diameters  large  compared  with  aggregate  sizes, 
and  also  widely  heterogeneous,  >  2,  distributions  do  not  conform  to 
those  of  the  model.  Also,  the  distributions  predicted  by  the  model  will 
not  apply  to  aerosols  that  have  undergone  significant  amounts  of  coagula¬ 
tion,  sedimentation,  or  diffusion,  after  the  initial  processes  of  genera¬ 
tion  are  complete.  Further,  if  in  the  generation  or  grinding  process, 
probabilities  for  interaction  with  breakup  stress  are  weighted  in  favor 
of  certain  size  fractions,  the  basic  premise  of  equal  access  to  the  stress 
by  all  particles  is  contradicted,  and  the  model  would  not  apply.  For 


ze  Distribution  Freeze-Dried  S.  marcescens .  Ground 
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example,  if  is  the  probability  that  a  particle  subjected  to  a  constant 
stress  breaks  up  in  a  fixed  time  interval,  tj_,  then  the  probability  that 
it  will  break  up  in  E^  At^  =  t^  is  given  by 


Absolute  size  distribution  predictions,  based  on  particle  properties 
and  dynamic  properties  of  the  compressible  flow  disseminating  system,  were 
successful  within  the  accuracy  of  model  for  the  one  powder  tested  and  on  a 
single  scale  of  dissemination.  In  this  study,  there  is  indication  that 
a  generally  applicable  number,  K"^-,  may  apply  for  aerodynamic  breakup  of 
solids,  similar  to  the  critical  Weber  number  for  liquid  drop  instability. 

An  upper  bound  on  the  number  was  obtained.  More  research  on  basic  particle 
contact  forces  and  bulk  tensile  strengths  of  powders  will  be  required  to 
establish  the  validity  of  this  concept,  however.  It  may  also  be  found  that 
K-1  is  a  function  of  the  compressible  flow  system,  involving  the  length 
of  time  an  aggregate  of  basic  particles  is  subjected  to  the  deagglomerating 
stress. 


For  friable  solids,  characterized  by  a  quasi-crystalline  structure, 
low  shear  stress,  and  high  elastic  modulus,  grinding  data  have  been  obtained 
supporting  the  concept  that  grinding  probabilities  are  proportional  to  the 
square  of  the  particle  diameter.  In  no  sense  can  this  be  applied  univer¬ 
sally  to  all  solids,  however.  For  example,  an  increase  in  size  of  particles 
of  Kaolinite  (refined  Devon  China  clay)  has  been  observed  during  ball 
milling.23’ 23  Further,  white  blasting  sand  demonstrated  an  increase  in 
particle  size  as  a  result  of  dry  grinding  by  mortar  and  pestle.  For 
materials  of  these  types,  at  least,  a  suitable  grinding  model  would  have 
to  predict  fusion  of  basic  particles  instead  of  fracture  at  some  point 
in  the  grinding  process. 
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